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We study the dynamics of a Heisenberg-XY spin chain with an unknown state coded into one 
qubit or a pair of entangled qubits, with the rest of the spins being in a polarized state. The time 
evolution involves magnon excitations, and through them the entanglement is transported across 
the channel. For a large number of qubits, explicit formulae for the concurrences, measures for 
two-qubit entanglements, and the fidelity for recovering the state some distance away are calculated 
as functions of time. Initial states with an entangled pair of qubits show better fidelity, which 
takes its first maximum value at earlier times, compared to initial states with no entangled pair. In 
particular initial states with a pair of qubits in an unknown state a || +/3 H are best suited for 
quantum state transport. 



Quantum entanglement has been recognized as 
an important resource for quantum information and 
computation^, for transmission of quantum states 
through a channelQ- There have been many propos- 
als of physical systems to serve as channels for quantum 
communication 3] , and in particular the spin chains 4]. 
The essential idea is to encode one particular qubit (spin- 
1/2 degree of freedom), and let it be transported across 
the chain to recover the code from another qubit some 
distance away. The quantum spin chains are well suited 
for state transport, as one could use the Schroeding dy- 
namics to propagate entanglement, and thus achieve the 
desired quantum communication. 

The dynamics of entanglement of spin systems, viz. 
the study of time evolution of an initial state, can be 
classified into three categories. The dynamics of a 
definite state or a S'^ non-definite state under con- 
serving time evolution (for example the Heisenberg-XY 
model) , and the time evolution of an initial state under a 
nonconserving dynamics (for example the transverse 
Ising model). The structure of entanglement sharing, the 
time scales for entanglement transport are very different 
between these categories. In this letter, we will study the 
Heisenberg-XY model, where is conserved through the 
time evolution. The main features of the dynamics of en- 
tanglement that will be addressed are, the appearance of 
pairwise entanglements between distant spins, viz. con- 
currences, starting from an initial state with no entan- 
gled pairs or exactly one entangled pair of spins, and the 
time scales for the transport of entanglement. The max- 
imally entangled initial states ti i 11 (the Bell state 
Bl) and TT ± ii (the Bell state B2) evolve quite differ- 
ently as we shall see below. The dynamics of entangle- 
ment in these states involve one-magnon and two-magnon 
excitations. The entanglement sharing in one-magnon 
and two-magnon eimstates of the Heisenberg-XY model 
has been studied0,Q. The time evolution of Bl states 
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involves the one-magnon excitations, which are not af- 
fected by — interactions, where as for the B2 states 
two-magnon excitations are involved, which introduce 
many complications. Here, not only two-magnon scat- 
tering states (equivalent to two one-magnon excitations) 
that have a weak signature of the interactions, but two- 
magnon bound states with a strong interaction effects^ 
have a significant contribution for entanglement dynam- 
ics. 

Let us consider an anisotropic Heisenberg-XY model 
for a linear chain of spins (s=l/2), with a Hamiltonian 

H = K,Y. - f E +h.c-by,si-e^ 

(1) 

where Kz,K are the interaction strengths for the z- 
components and the xy components respectively, and B 
is the strength of the magnetic field along the z direction. 
The constant = -NKj4,-NB/2 is the energy of the 
ferromagnetic state, with all the spins poplarized along 
the z direction. The above Hamiltonian generates uni- 
tary evolution from the Schroedinger equation which con- 
serves the total S*^ in the state. The energy eigenstates 
and eigenvalues are known exactly through the Bethe 
ansatzj^. The ferromagnetic state \F) = \] .. ]) has no 
dynamics, being an eigenstate with zero eigenvalue, and 
has no entanglement between any pair of spins. Let us 
first consider an unentangled state with one spin, at the 
site I, in an unknown state, at time t=0 given by 

= I T .. T>|a T +/5 i)i = {a + Psi)\F). (2) 

This state has no entanglement, being a direct product 
state of different site states. Let us denote \n) = s~|i^), 
a state with one down spin at site n. At a later time t, 
the state can be written as 

|^„(<)) =a|F)+/3^0„(t)|n) (3) 

where 0„ = \/N'Y^^e'x:pi{q{n — I) — Eqi/h), and the 
one-magnon energy Eq = —Kcosq + B. In the limit 
of large number of spins, N oo, the wave function can 
be expressed in terms of the Bessel function J„_; as 

Mt) = e-^^e'^("-') J„_z(ifi/^)- (4) 
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The time scale for the structure in the wave function is 
T = h/ K , and from now on we will write the time as a 
multiple of r, as T = i/r. We can estimate the time scale, 
using K ^ O.Oley, as t ~ 10"^'^ sec. And the magnetic 
field adds on a constant phase to the wave function, and 
thus can be dropped. At time t the mixed state of a 
given site i can be denoted by the reduced density matrix 
Pi = tr[ ^|V'«)('0u|, where the prime indicates a partial 
trace over all states except at site i. It is straightforward 
to write down the reduced density matrix as 



= 



'\Mm\r){u + w\'\Mt)ni){i\- (5) 



Now the initial unknown state encoded in th I'th qubit, 
can be extracted from the i'th qubit with a fidelity Fi — 
TrpuPi where p„ = |a T +P T +P I \, which works 
out to be (with i — l+r, a distance r away from the initial 
site)F,(T) = |a|2 + |/?|2(|/5|2-|a- 
over all possible initial states, i.e. 
average fidelity is 



^)J^(r). By averaging 
the Bloch sphere, the 



(6) 



A similar formula has been derived in (J| , except our cal- 
culation is simplified due to the limit of large N. Now, we 
are interested in propagating the initial code to a distance 
r, and recover it. The fidelity Fr has a maximum value 
for T K r, which means the quantum state is transported 
at a rate Vt = l/r. After waiting for a time t = Vtr we 
have the best recovery of the quantum state at site a dis- 
tance r away from the initial site /. In Fig.l, the fidelity 
has been plotted as a function of time, for r — 100; the 
first maximum is for Tc ~ r, along with the result for 
entangled initial states we shall discuss below. 

Though the initial state has no pairwise entanglement, 
for T 0, the state develops entanglement. We use the 
concurrence measure 7] for the pairwise entanglement, 
which can be calculated from the two-site reduced den- 
sity matrix pij (which is obtained by tracing over all spins 
except those at sites i and j). The time-reversed density 
matrix is denoted by pij, and the eigenvalues of pijPij by 
A1..A4 in the descending order. Then the concurrence be- 

f3 - - - - 

tween the two sites is,7j, Cy = max(Ai^ — A| — A| — A4 , 0). 
Here, the two-site density matrix has the form 



/I 



V 



I fl|2 



given by [El C, 



0/ 



Now the concurrence 

2\l3\'^\Ji_i{T)Jj_i{T)\. A plot of concurrences vs 
T are shown in Fig. 2, for i = I + l,j = I and 
i = I + 2,j = I. For small T, the concurrences grow 
as Cij w= {2\(3\'^/r\){T/2Y. For large T, we have 

di w with oscillations. 

Now we will turn to initial states with entangled pairs 
of spins. Entangled states are expected to better than 
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FIG. 1: The average fidelity Fr as a function of T, for r = 100, 
for the three types of states discussed. For the entangled 
states, s — I — m — 25 has been used. The first maximum for 
the unentangled is at T = r = 100, and for the Bl and B2 
states, at T = r — s. 



the unentangled states considered above. The dynamics 
will transport and further generate entanglement, as we 
shall see below. Let us first consider an initial state with 
a pair of qubits at sites I and m in an entangled state (Bl 
state) a\ ti) -I- P\ it), all other spins polarized, which is 
represented as 

|Vi(r = 0)) = (asr+/3s-J|^^) =^0„(O)|n). (7) 

For T = all concurrences are zero except, C;„i — 1. For 
T 7^ 0, this entanglement spreads, and is transported to 
other pairs. Again using the one-magnon excited states, 
we can write down the wave function as a function of T 

0„(T) = /3e't("-') J,_z(T) + ae^^f"-") J._™(T). (8) 

Since the above state is a one-magnon state (though may 
not be an eigenstate) the concurrence is given by|5| Cij = 
2\(j)*(j)j\. In particular, the concurrence between the sites I 
and m for later times, for the maximally-entangled initial 
state {a = P = 1/ V2) , for ^ — TO even or odd respectively, 
is given as 

G™ = (Jo-f (-l)'^Ji-™)^ or {J^ + JlJ. (9) 

In Figure 2, C;™ has been plotted as a function of T, for 
I — m ^ 1 , which drops from the initial value of unity as 

Now, the fidelity of recovering a quantum state \a t 
+f3 i) at a site i can be calculated straightforwardly, 
analogous to the unentangled state we calculated before, 
as (for i — m + r, a distance r away from the site to) 
F^iT) = lap + (|/3p - \a\^Mr+m{T)\^. For \a\ « \f}\, the 
fidelity is close to 1 /4, as before for the unentangled case, 
which means the channel is noisy for a good recovery. In 
this case the unknown state can be recovered from two 
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FIG. 2: The concurrence dj as a function of T, using — 
1/2. It can be seen how concurrence builds up between the 
initial site I, I + 1 and I + r for the unentangled initial state, 
and how the concurrence decreases between the initial sites I 
and m for Bl and B2 Bell states. 

sites, as we shall discuss below. The average fidelity, after 
averaging over the Bloch sphere, is given by 

Fr{T) = \ + \{Jl_^i_„,^~Jl) (10) 

which should be compared with the formula we obtained 
(Eq.6) for the unentangled initial state. The first maxi- 
mum of the fidelity now depends on Jr-(i-m), which will 
occur for an earlier time, and the maximum value will be 
more than for J^. The rate at which the state propagates 
is still vt ~ 1/t, but we need to wait for a shorter time 
interval t ~ vt{r ~ (l — m)) to recover the state a distance 
r away from the m'th site. The distance between the 
two initial sites can be chosen conveniently for a given 
distance over which the transport is desired. The fidelity 
for this case has been shown in Figure 1 for a few values 
of r as a function of time. Due to the presence of the 
competing term (the last term in the above), the fidelity 
could fall below the value for the unentangled state, after 
the first maximum. 

The initial entangled state at sites I and m,pBi = \a ti 
+f3 11) {a ti +/3 it): can be extracted from sites i and j 
at a later time with a fidelity, Gy = TrpsiPij , which can 
be calculated after some manipulations simply as Gy — 
\a(f)j + (3(j)i\'^. This has a maxima structure even for |a| ~ 
unlike the fidelity function Fi discussed earlier. For 
i = l + r,j — m + r, i.e. for a pair of sites translated 
by r from the initial pair l,m, the fidelity takes the form 
(using s ~ I — m), Gr — \2al3Jr + o? exp {iirs /2)Jr+s + 

exp (— i7rs/2) Jr-sp. After averaging over the Bloch 
sphere, it has a simple form 

In Fig. 3, we show the fidelity Gr as a function of T for 
r = 50,100. 



FIG. 3: The average fidelity Gr as a function of T, for r = 
50, 100 using s = Z — m = 10. A constant 0.5 has been added 
on for Bl state, to show both on the same graph. 

Let us now turn to the most difficult case of an initial 
state with a B2 state at sites / and m, and the rest of 
the spins polarized, given as {a + /3s;~s~)|F) = a\F) + 
f3\^{T = 0)). The state can be written as 

\MT = 0)) = a\F) + pJ2^,jiO)\ij), (12) 

where \ij) stands for a two-magnon basis state with two 
down spins at sites i and j. Initially for i — I, j = m the 
wave function is unity, and zero for all other values of i 
and j. The time evolution of the second term above can 
be worked out in terms of the two-magnon excitations, 
using the Bethe Ansatz solution for two down spins. The 
concurrences in the two-magnon eigenstates have been 
worked out 6] both for the scattering and the bound 
states. The wave function is extremely complicated due 
to the magnon interactions arising for a nonzero Kz, 
which can be vastly simplified by taking the limit K))Kz, 
i.e. dropping the interaction terms in the Hamilto- 
nian. In the XY limit, the wave function as a function of 
time takes the form (after taking N ^ oo limit) 

4j(T) = e^^(*+^-'-™)(J,_,(r)J,_,„ - J^-,nJj-i), (13) 

which is antisymmetric in the two indices, reflecting the 
underlying fermionic nature of the moving down spins. 
Now, following through the steps as before, the two-site 
reduced density matrix has the form 

/ |a|2+/3|2u,, a(3*cl>*A 

In the above, the various matrix elements stand for Uij = 
((l/2 + sf)(l/2 + s|)),z;,, = ((l/2-sf)(l/2-s|)),z«H, = 
((l/2-sf)(l/2 + s|)),u;2„- = ((l/2-f<)(l/2-s|)),z,, = 
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(s^s~), where the expectation value is taken in the two- 
magnon state \^{T)) only. Now following through the 
further steps of Wootters0, the eigenvalues of PijPij 
for the above density matrix are (^/(laP + ± 
\a* ,\f3\'^{y/wiW2 ± This gives two regimes 

for the concurrence as 

= 2\a*f3(j)^j I - 2\P\^^wi,jW2rj (14) 

which ever term is positive, and otherwise zero. The off- 
diagonal matrix element can be calculated as (taking 

J + 1 

where the sum in r]ij is over all values of n, which 
can be calculated using the addition rule Jn{x + y) = 
X]fc Jk{x)Jn-k{y)i Eind settig the site m at the middle of 
the chain for convenience, 

r]ij = e'^^J-'^ [Jr-lJj-l + J^-mJJ-m)■ (16) 

And , which is just the finite sum, is quite complicated 
to calculate in general. The diagonal matrix elements 
are = ,Uij « l,wiy = rju - |(/)ijP,W2ii = Vjj - 
This simplifies the expression for the concurrence 
between two sites i and j as 

Cy- = |/3|max(O,|/3||z|-|0|,2(|a0|-|/3|>^)). (17) 

For I — m odd, the off-diagonal matrix element zi„i — 0, 
and the expression for the maximally- entangled initial 
state is simple, Cim — Jq + Jf-jm which is what we got 
for Bl Bell states. For / — m even, the expression is 
still complicated. For I — m ~ 2, Cim — max(0, \z\/2 — 
\cj)\/V2, |(/)| - wi), where = (J^ - j|), = + j| - 
\(j)\'^,z = 2J0J2 + Ji (Jo + >/2)^- The concurrence between 
the sites I and m is plotted as a function of T, for l—m=l 
in Fig. 2, along with the result for the Bl states and the 
unentangled state. 



The fidelity of recovering the state \a t +/3 i) at a site 
i is again straightforward, Fi = \a\^ + |/3p(|/?p — lap)?]^^. 
The average fidelity, for i = m + r, 

Fr = \ + \{JUl-rn)+Jr)- (18) 

As compared to the expression for the Bl states, in the 
above there is no competition between and Jr-(i-m)- 
The fidelity here is greater than that of the unentangled 
state for all times. The first maximum value is deter- 
mined by the first term, as in the case of Bl states, which 
occurs ioY t — vt{r — I -\- m). A comparison of the fidelity 
as a function of T for all the three cases, from Fig.l for 
r = 100, shows that the B2 states have better fidelity. 
Now, the fidelity of recovering \a IT +/5 ii) from sites i 
and j, analogous to the Bl states we discussed before, is 
given as G,, = \a\^ + + \a(}\^{r^j + '/>.j)- This 

function also exhibits a maxima structure for |a| « \(3\. 
The average fidelity is (for i = I + r, j = ni + r) 

Gr = \ + \{j'^ - Jr-sJr+s){Jr ~ Jr-sJr+s + e'^^). (19) 

The fidelity as a function of T is shown in Fig. 3, for 
r — 50, 100, s = 10, along with the result for Bl states. 
The B2 states exhibit better fidelity here also, as is the 
case in Fig.l. 

In conclusion, we have investigated the quantum state 
transport across a channel of qubits, the spin chain, using 
the Heisenberg-XY dynamics. The presence of entangle- 
ment and its dynamics is crucial for communication over 
the channel. Initial states with a pair of qubits in a state 
a 1 1 +(3 [[ show better fidelity. Here, it will be interest- 
ing to investigate the effect of a nonzero K^', significant 
changes in the wave functions, the concurrences and the 
fidelity are expected. Finally, states with many entan- 
gled pairs in an optimized network may demonstrate an 
almost-ideal quantum communication, that is, a telepor- 
tation protocol with the sender and the receiver at a fixed 
distance and the code transported with negligible inter- 
ference from the network channel. 
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